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Abstract. The dynamical evolution of the Gould Belt has been modelled in 3D and confronted to the spatial 
and velocity distributions of all HI and H2 clouds found within a few hundred parsecs from the Sun and to the 
Hipparcos distances of the nearby OB associations. The model describes the expansion of a shock wave that 
sweeps momentum from the ambient medium. It includes the effects of the Galactic differential rotation and its 
gravitational torque, as well as interstellar density gradients within and away from the Galactic plane, possible 
fragmentation and drag forces in the late stages, and an initial rotation of the system. The evolved Belt geometry 
and velocity field have been fitted to the directions and velocities of the nearby clouds using a maximum-likelihood 
test. In order to do so, local clouds have been systematically searched for in the available HI and CO surveys. 
The likelihood function also included a distance estimate for a subset of well-known clouds. 

The best fit to the data yields values for the current Belt semi-axes of (373 ± 5) pc and (233 ± 5) pc, and an 
inclination of 17.2° ±0.5°. These characteristics are consistent with earlier results, but a different Belt orientation 
has been found because of the presence of new molecular clouds and the revised distance information: the Belt 
centre currently lies (104 ± 4) pc away from the Sun, towards the Galactic longitude Icmtre ~ 180.4° ± 2.0°, and 
the ascending node longitude is In = 296.1° ±2.0°. The Belt characteristics are independent of an initial rotation. 
The present Belt rim is found to coincide with most of the nearby OB associations and H2 clouds, but the Belt 
expansion bears little relation to the average association velocities and the younger ones are surprisingly found 
farther out from the Belt center. An initial kinetic energy of (1.0±0.1) 10*^ J and an expansion age of (26.4±0.4) 
Myr are required, in good agreement with earlier 2D estimates. The factor of 2 discrepancy that exists between 
the dynamical Belt age and that derived from photometric stellar ages could not be solved by adding a vertical 
dimension in the expansion, nor by adding drag forces and fragmentation, nor by introducing an initial rotation. 
Allowing the Belt to cross the Galactic disc before reaching its present position would require a longer age of 52 
Myr, but the very poor fit to the data does not support this possibility. 

Key words, shock waves - ISM: clouds - ISM: kinematics and dynamics - open clusters and associations: individual 
(Gould Belt) - solar neighbourhood 



1. Introduction 

The Gould Belt is a nearby starburst region where many 
stars have formed over 30 to 40 million years in a sur- 
prisingly flat and inclined disc. New facets of its activity 
have recently emerged at high energy with the dis covery of 
a population of 7-ray sou rces associated with it llGrenieil 
l200nt ICxehrels et al.l l200rfl. No clear picture has emerged 
yet as to the nature of these objects. Neutron star ac- 
tivity in various forms appears as a promising prospect, 
in particular 7-ray emission from million- year old pulsars. 
As supernova relics, these sources have drawn attention to 
the enhanced supernova rate inside the Gould Belt, which 
has been estimated to be 3 to 5 times higher that in the 
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local Galactic disc. In other words 20 to 27 supernovae 
have occurred in the Belt per m illion year over the past 
fe w millio n years (Grenier"200(/). 

iHerschel (1847 ) first pointed out that the distribution 
of bright stars was asymmetric about the Galactic plane. 
The geometry of this young structure (most of the stars 
ar e less th an 30-60 million years old) was studied in Il874l 
bv lGouldl who determined its inclin ation and the d irection 
of its poles. We refer the reader to IPoppell l)l997j) for an 
extensive review of Gould's Belt system and its relation to 
the local interstellar medium. Besides its peculiar geom- 
etry, this system of early-type stars is known to expand, 
and to rotate in the same direction as the galactic rota- 
tion. The Gould Belt also contains interstellar clouds and 
iLindbladI lll967ll gave strong evidence for the Belt relation 
to a locally expanding HI ring. Famous H2 complexes, such 



2 



Perrot,C.A. and Grenier I.A.: 3D dynamical evolution of the Gould Belt 



as Orion and Ophiuchus, are often mentioned in relation 
to the Belt and the fact that nearby dark cl ouds partici- 
pate t o the Belt expansion was recognized bv lTavlor et al.1 
l|1987(l . Since then, the molecular components of the lo- 
cal interstellar medium have been thoroughly mapped in 
the CO surveys, but no correlation study has been under- 
ta ken. 

IComeron et"al] l|l994|) concluded that 40% to 50% of 
the young massive stars (with spectral types earlier than 
B8) that lie within 450 pc from the Sun belong to the 
Belt. The new Hipparcos data bring this fraction to 60% 
for stars within 60 pc and a Belt age of 30 to 60 Myr 
l|Torra et al.ll2000|) . Whereas the massive-star content of 
the Belt has been extensively studied, less is known about 
its low-mass star production. Young (30-80 Myr old) 
Lithium-rich solar-mass stars with active coronae show 
up as X-ray sources and nicely trace the Belt in the sky 
i Gtiillout et alJll998(l . Even though the X-ray horizon is 
limited by interstellar absorption to 100 or 150 pc for these 
faint sources, their space distribution indicates that stellar 
formation is not only active along the Belt rim, but also 
inward, over a significant, yet poorly constrained, radial 
extent. 

Preserving the structural coherence (flatness and incli- 
nation) of the Belt stellar system over a long time span is 
challenging, especially over a large fr action of the ver tical 
oscillation and expansion timescales. IComeronI |l999) has 
explored the stellar vertical motions and found that the 
rotating stellar disc was initially tilted. Its rotation axis 
was not perpendicular to the Galactic plane as would be 
expected from the disruption of a giant ro tating molecular 
cloud. According to lLindblad et al.l 1 19971) . the dissolution 
of an unbound rotating system of stars, possibly born 30 
Myr ago in a spiral arm, may reproduce the observations, 
and the rotation may explain the persistence of a flat disc. 

The Belt flatness and its tilt, some 20° to the Galactic 
plane, bear important information on the Belt origin, but 
remain very difficult to interpret. Various scenarii involve 
the oblique im pact of a high-ve locity HI cloud on the 
Galactic disc (,Comer6n fc Torral 11992. .19941') or acas- 
cade of several supernova explosions (see 
for a review). The measured Oort constants of the stel- 
lar field are consistent with a cl oud impact about 50 
Myr ago ijComeron fc Torral 1 1994j) . On the other hand, 
an explosive event would explain the kinematics of the 
cold neutral medium at |6| > 10° ( Poppel fc Marronettil 
studied the expansion, within the 
Galactic plane, of an initially circular shock wave. The 
30 Myr expansion in this 2D model was constrained by 
the present-day longitude-velocity distribution of the HI 
gas in the Lindblad ring. Following the expansion of a 
supe rbubble, with or wit hout ambient interstellar pres- 
sure, iMoreno et alJ l|l999t) compared the emerging stellar 
orbits to the velocity field of nearby massive stars. They 
concluded that the Pleiades group, and possibly the Sco- 
Cen association, significantly d eviate f rom t his single ex- 
pansion model. More recently, lOlanol l|200ll) presented a 
3D model for the local system of gas and stars associ- 



ated with the Sirius supercluster, the Gould Belt, and the 
Local Arm. The latter two subsystems would have formed 
in the braking of a supercloud in a spiral arm while the 
stars of older generations, as in the Sirius supercluster, 
would move on with their initial velocity field. The su- 
percloud angular momentum being concentrated at large 
radii, the inner regions would collapse into a flattened disc, 
precursor of the Gould Belt, whereas the ejection of the 
outer parts into a super-ring would form a precursor of 
the Local Arm. 

To improve on Olano's model and to study the Belt re- 
lation to the new molecular clouds and its position with re- 
spect to the local clouds and OB associations, we model in 
this work the 3D dynamical evolution of an inclined cylin- 
drical shock wave, expanding in the interstellar medium in 
the Galactic plane and at higher altitudes (c./. section O. 
Using a maximum-likelihood approach described in sec- 
tion 12.41 we confront the characteristics of the gas shell 
in the longitude-latitude-velocity phase space to the ob- 
served positions and velocities of all nearby HI and H2 
clouds that have been found in extensive HI and CO sur- 
veys (see section|31l. We then discuss the best-fit evolution 
and current geometry of the Belt and its relation to the 
stars in section 0] 

2. Evolution scenario and analysis 

2.1. Dynamical evolution 

The model does not attempt to explain the origin of 
the outburst, but describes the lateral expansion of an 
inclined, cylindrical shock wave that sweeps momentum 
from the ambient interstellar medium. As a consequence 
of the Galactic differential rotation, the circular section of 
the Belt rapidly evolves into an elliptical one. Additionally, 
the combined actions of the Galactic gravitational poten- 
tial and of the interstellar density gradients further warp 
the Belt. The Belt ring was split into 60 elementary sec- 
tions. Their contiguous surfaces delineated the Belt rim 
without any holes at each time step and allowed it to 
bend and warp easily. The motion of each Belt element 
was simulated, taking into account the galactic gravita- 
tional forces, and the momentum variation induced by the 
swept-up gas. The integration was done using a Runge- 
Kutta algorithm with an adaptative timestep, so that the 
distance step for the fastest gas element is less than 1 pc. 

The initial Belt geometry is described by a limited 
set of free parameters: its inclination with respect to the 
Galactic plane, Lpf)\ the longitude of its ascending node, 
Zjj; the constant cylinder height, H\ the longitude and dis- 



tance from the Sun of its centre, dPcentre ^-^d ^centre^ 



its 



age, r; the initial mass and velocity of the ejecta, Mq and 
vq. The expelled mass was uniformly distributed along the 
Belt. An initial radius of 20 pc, typical of OB associations, 
was adopted, but it does not significantly influence the 
subsequent evolution. The Belt centre was assumed to be 
initially located in the Galactic plane, at 2; = 0. An ellip- 
tical fit to the position of the sections allows to describe 
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the Belt geometry at any stage in terms of its semi-major 
axis, a, and semi-minor axis, b. 

2.2. Interstellar density gradients 

The local HI and H2 gas density into which the shock wave 
expands is described by a com bination of three Gaussian 
and one exponent ial functions I Dickev fc LockmanlllQQCt 
lDa,me et a,1.llT987i): 



Table 1. Longitude, latitude and velocity ranges for cloud 
selection 



+N-i,Hi exp[ ) + exp 
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with Ni,Hi = 0.395 cm^^ di = 90.03 pc, N2,hi = 0.107 
cm-3, (72 = 225.17 pc, iVg,^/ = 0.064 cm'^,' erg = 403.0 
pc, = 0.2 cm~'^, and (7h-2 — 74.0 pc. A radial, galac- 
tocentric density gradient was also included, dropping lin- 
early by a factor of 2 over 1 kpc. It had no impact on 
the results, so the sole z dependence on the interstellar 
density has been retained hereinafter. In order to com- 
pute the interstellar gas momentum, we adopted the lAU 
recommended values for the Oort's constants Ac ~ 14.5 
km s~^ kpc~^ and Be = —12.5 km s~^ kpc~^, the Solar 
galactocentric radius Rc = 8.5 kpc, and its orbital velocity 
Vc = 220 km s~^. No internal pressure from supernova or 
stellar winds was added, nor any external pressure from 
the interstellar medium. All calculations were done in a 
Galactic cartesian inertial frame, centred on the Galactic 
centre, the x-axis pointing from the Sun to the Galactic 
centre ai t — Q, the y-axis pointing in the direction of 
Galactic rotation, and the z-axis pointing to the North 
Galactic pole. 
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compared with the space and velocity distributions of the 
nearby HI and H2 clouds. Because the direction of these 
clouds is known, but rarely their distance, we developed a 
specific maximum-likelihood test to compare their direc- 
tion and radial velocity with those of the modelled Belt 
sections. The likelihood expression is based on the prob- 
ability "^ij for a cloud i to be seen at relative distance 
Dij from an elementary Belt section j, and at a relative 
radial velocity Vij with respect to that of the Belt section 
(c./. Equation |3J| . Gaussian probability density functions 
have been adopted in relative distance and relative radial 
velocity with standard deviations <jd = 20 pc and Uy — 1 
km s~^, respectively. So, is given by 



2.3. The Galactic gravitational potential 

The torque locally induced by the gravitational poten- 
tial of the Galactic disc was calculated as a function of 
altitude z using the local stellar mass density distribu- 
tion, with a volume den sity = 7.6 10"^ pc~^ at 
z = ijCreze et all 19981) and an exp onential scale height 
= 260 ± 60 pc ijOiha et alJll996|) . Assuming that the 
Galactic disc can be locally described by an infinite plane, 
the corresponding vertical acceleration is given by: 



"Hi, ^ 



4 7r2 a3 



(3) 



D 



where 1/gs is the Belt perimeter. 

When the cloud distance to the Sun, Di ± ADi, is 
available (c./. Tables IXtI and IX!^ . the cloud likelihood is 
obtained by integrating over the 60 Belt sections: 



1 — cxp ( 



z 

N 



at = -47rC/p*z* 
Q denoting Newton's gravitational constant. 



(2) 



2.4. Maximum-likelihood analysis 

Varying the initial expansion parameters, the geometry 
and velocity field of the evolved Belt (today) have been 



J Belt 



Df Arj, AD, Aw,- dsi 



(4) 



Afti = Ali Ahi cos&i represents the 1 a confidence solid 
angle around the direction [li, bi) of peak emission in the 
cloud, in Galactic coordinates. Avi notes the 1 a error on 
the observed cloud centroid velocity. 

For a cloud with no distance estimate, we performed a 
Gaufi-Legendre integration of the spatial probability den- 
sity function along the cloud {U, bi) direction: 
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Fig. 1. Galactic longitude and latitude distribution of the nearby HI (black dots) and H2 (white dots) clouds. The 
dots mark the direction of peak emission. Clouds at \b\ < 5° or 8° have not been selected. The curves follow the trace 
of the Belt across the sky as obtained from the best fit. The upper and lower lines mark the corresponding edges of 
the 60 pc thick disc of the Belt. 
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Xj , Uj , and Zj note the cartesian coordinates of the centre 
of a Belt section. 

The likelihood C{(po,l^, H, d", r, Mq, vq) of a model 
is obtained by the product of individual (indepen- 
dent) cloud probabilities Ci. The maximum likelihood is 
found using the down hill simplex method developed by 
iNelder Me^ lllflfi,^. The iteration stops when the dif- 
ference between the function values at each point of the 
simplex and the current minimum are less than 10"'\ This 
procedure has been successfully tested using simulated 
clouds randomly generated from a known Belt. 

3. Cloud selection 

As noted by several authors, the Gould Belt has left its 
imprint on the local gas dynamics as seen in the HI data. 
The local molecular clouds, as mapped in CO, appear to 
follow a longitude-velocity patter n similar to that o f the 
atomic gas in the Lindblad ring (iDame et a,lJll98'i^ . To 
include all nearby atomic and molecular clouds, we used 



th e extensive CO sur vey of the Milky Way at |5| < 25° 
bvlDame et all lll987l) and two complementary H I surveys 
i Hartmann fc BurtonI Il997t IStrong et alJll982|) to cover 
the whole sky. Clouds were systematically searched for 
and isol a.ted using the p owerf ul CLUMPFIND tool devel- 
oped by IWilliams et al This algorithm searches 
for coherent peaks of emission based on closed contours in 
a longitude-latitude-velocity data cube, from the highest 
intensity to the lowest. The algorithm yields the position 
of peak emission and centroid radial velocity of the clouds. 

The CO data cube was processed in the [—15, +17] 
km s~^ velocity range characteristic of the local medium. 
The temperature increment for the CLUMPFIND to fol- 
low down the intensity levels was set to AT = IK. 
The same procedure was applied to the HI data cubes 
for velocities |u| < 20 km s~^ and latitudes restricted 
to |6| > 3° to leave out extensions to medium latitudes 
of background clouds in the Galactic plane. Temperature 
steps o f AT = 0.07 K and AT = 2 K ha ve been used 
for the IHartmann fc BurtoiJ l)l997() and the IStrong et alJ 
(ill98^ surveys, respectively. 

Table n summarizes the selection criteria used in lon- 
gitude, latitude and velocity to avoid contamination from 
background clouds in the Galactic disc or in the Local 
Arm. These cuts were kept minimal to limit biases and 
to avoid an a priori selection of clouds obviously linked to 
the Belt while retaining the intrinsic spatial and velocity 
dispersion of the clouds in our neighbourhood. The lists 
of selected HI and H2 clouds and their characteristics are 
given in Tables FOI and 1X^21 respectively. The (l,b) direc- 
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Fig. 2. longitude- velocity distribution of the nearby HI 
(black dots) and H2 (white dots) clouds and the best fit 
to the Belt rim kinematics. Projected velocities are given 
in the Local Standard of Rest. 
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Fig. 3. longitude, latitude, and velocity distribution of the 
nearby HI (black dots) and H2 (white dots) clouds and the 
best fit to the Belt rim kinematics. Projected velocities are 
given in the Local Standard of Rest and are colour-coded 
from black to white from negative to positive values 



tions and velocities of the selected clouds are displayed in 
Figures ^ and 13 



4. Results and discussion 

The initial and current characteristics of the Belt expan- 
sion that best fit the present-day cloud data are given 
in Table El The resulting geometry and dynamics appear 
to be independent of the presence or not in the fit of the 



Taurus clouds which, since they lie near the centre, do not 
participate to the expansion of the Belt rim . The quoted 
la errors are purely statistical and have been determined 
from the log- likelihood ratios and the information matrix 
llStronelllflsi . 

Figure ^ shows the maximum-likelihood evolution of 
the Belt rim, as seen at different epochs in a plane per- 
pendicular to the Galactic one. The current Belt geometry 
is close to that depicted in the T = 30 Myr plot and the 
last plot illustrates how the Belt may fall back onto the 
Galactic disc within 10 to 15 Myr. The height that best 
fits the data is = 60 ± 1 pc. The expansion proceeds 
faster at higher z altitudes because of the reduced ambient 
interstellar density. The inclined Belt starts by rapidly ex- 
panding out of the Galactic disc. Its shape gradually nar- 
rows into an ellipse that precesses because of the swept-up 
momentum from the interstellar gas and its differential ro- 
tation. The inclination also decreases with time and the 
disc gets clearly warped because of the gravitational pull 
from the Galactic disc. 

To zeroth order, the expansion radius R(t) is well de- 
scribed by momentum conservation in a uniform medium 
of density p: -^{npHR^R) — 0. The resulting power- law 
dependance R oc t^^^ is found in very good agreement with 
the evolution of the average size, R(t)=(a-|-b)/2, obtained 
from the model after the first 200 kyr. 

The best fit yields a position and velocity field for the 
present Belt that are in good agreement with the cloud 
data, as illustrated in Figures and 13 Figure ^ shows 
the trace of the Belt across the sky in longitude and lat- 
itude. Apart from a few points, which have not been re- 
moved from the sample to avoid biases, the modelled Belt 
geometry nicely fits the data. Its angular width increases 
at low longitudes where the rim is closer and the vary- 
ing angular width matches the angular dispersion seen 
in the nearby clouds. Figure |21 shows the cloud distri- 
bution in longitude and velocity, as well as the velocity 
field of the present rim. The slow expansion of the Belt to 
date reproduces the double-peaked velocity pattern over 
all longitudes. The velocity pattern clearly differs from the 
vlsr = A.d.sin{2l) .cos^ (b) dependence expected from a 
pure differential rotation. The difference reflects the tilt of 
the Belt and its expansion. The peak-to-peak velocity am- 
plitude is, however, not quite satisfactory, even in the face 
of the large intrinsic dispersion in the cloud data. This 
aspect will be addressed in section 14.21 Figure El shows 
that the model reasonably reproduces the complex data 
distribution in the longitude, latitude, and velocity phase 
space. 



4.1. Geometry 

The present ellipse has a semi-major axis a ~ 354 ± 5 pc 
and a semi-minor axis b — 232 ± 5 pc, in good agree- 
ment with the dimensions of 360 by 210 pc and 341 by 



ment witn tne ai mensions 01 d bU b y ziu pc ana ^ 341 p 
267 pc derived by'Olano (1982^ and' iMoreno et all l|l99£ 
from the sole HI data when modelling the expansion of 
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Fig. 4. The Gould Belt evolution as seen at different epochs after the outburst, in a plane perpendicular to the Galactic 
plane, centred on the Belt centre. The x axis points to the Galactic centre and the location of the Sun is marked by 
an asterisk. The Belt expands and is further distorted by the combined effects of the Galactic differential rotation, 
interstellar density gradients above the Galactic plane, and the gravitational pull of the Galactic disc. The plot for 30 
Myr is representative of the present-day geometry. 



a 2D ring in the Galactic plane or a 3D superbubble. 
The maximum-likelihood fit also yields a present inclina- 
tion = 17.2° ± 0.3° on the Galactic plane that nicely 
compares with the stellar estimates. Inclinations of 22.3°, 
16° — 22°, and 20° have b een found fro m the ma ssive-star 
pop ulation bvlGomeron et al. (19941. .Torra et al.. (JOOO), 
and lOlancl I 200lh . resoectivelv. A lower value of 12.5° was 



derived from the molecular clouds l|Tavlor et al.lll98^ . 

The difference with the present result can be attributed to 
a different cloud dataset. Taylor et al. did not select high- 
latitude dark clouds that are shown below to be consistent 
with the Belt shell. An inclination near 20° supports a 
strong connection between the Gould Belt and the Local 
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Fig. 5. Present position of the Gould Belt projected on the Galactic plane. The x and y axes point to the Galactic 
Centre and in the direction of the Galactic rotation, respectively. Both are centred on the Belt centre. The velocity 
field outlines the Belt expansion with respect to the Local Standard of Rest. Nearby OB ass ociations are plotted 
as thick circles using Hipparcos estimates of their distance and dimensions from lde Zeeuw et al.. 1,1999) . The shaded 
circles mark the location of the main nearby H2 clou d comp l exes. The thick and thin ellipses note the Belt rim as 
obtained in this work and earlier from the HI data bv lOland l)l982l) . The triangle and star note the Belt centre and 
the Sun, respectively. 
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Fig. 6. 3D view of the present Belt wave and its velocity field ( with respect to the LS R) amid the local OB associations 
that are marked as spheres with a radius equal to their size l|de Zeeuw et al.lll999|) . The diamond and star note the 
Belt centre and the Sun, respectively. 
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Initial and present characteristics of the Gould 
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Bubble (or Chimney) since the axis of the latter appea rs 
to be perpendicular to the Belt plane l|Sfeir et alJll999|) . 

We find a position and orientation that significantly 
differ from previous estimates. The Belt centre is found 
at a distance rf*°^tre — ± 4 pc from the Sun, in the 
direction l*°^tre = 180° ± 2°. More distant centres, well 
int o the se cond quadr ant, were propo sed from the HI data 
by lOlanol ((^§2) and iMoreno et all (|l999) (131 pc and 
173 pc away, towards 166° and 117°, resp.). The ellipse 
axes also point to a different direction. Figure El gives the 
current Belt position projected onto the Galactic plane. 
The change in orientation can be attributed to the refined 
distance information used here and to the presence of new 
major nearb y complexes, suc h as Cepheus, Cassiopeia, 
and Polaris. iTavlor et al did not select their dark 

clouds as Belt members on a direction basis, but their 
combined direction, velocity, and distance appear here to 
be quite consistent with their being part of the expanding 
shell in the second quadrant (see Figures and |2l . 

Finally, the best fit yields an ascending node longitude 
In = 296.1° ± 2.0° that is - 10° higher than the values 
obtai ned from the stars { lo, = 284.5° from the massive 
ones llComeron et al.lll994l) and = 282° from the young 
low-mass stars ijCuillout et al.lll998l) ). The difference pos- 
sibly reflects the time-lag between the stellar birth epoch 
and the slowly precessing rim today. 

As shown in Figure the position of the shell nicely 
coincides with most of the nearby OB associations such 
as Per 0B2, Ori OBla & OBlc, LCC, UCL, US, Cep 
0B6, Tr 10, Vela 0B2, the po sitions of which are k nown 
from Hipparcos measurements l|de Zeeuw et al ■Il999l) . Col 
121 and Cep 0B2 clearly lie outside the Belt. Projections 
can be misleading. The 3D view displayed in Figure 
illustrates that Lac OBI, often quoted in relation to the 
Belt, lies on the wrong side of the Galactic plane to be part 
of it, even when taking into account the Belt thickness of 
60 PC . The old age of the as sociation, of the order of 12-16 
Myr ((de Zeeuw et alJll999l) . also confirms that Lac OBI 
is not related to the Belt. 



4.2. Effects of an initial rotation, stiell fragmentation 
and Galactic plane crossing 

The velocity dispersion of the data points in Figure |21 is 
large, but consistent with the cloud to cloud dispersion ob- 
served in the interstellar medium. A larger peak-to-peak 
amplitude in the shell velocity would, however, better rep- 
resent the data. The influence of the initial conditions, the 
Belt thickness, the distance information, and the interstel- 
lar density profiles have been investigated to reduce this 
discrepancy, but none of these options leads to a signifi- 
cant increase in velocity amplitude within sensible limits. 
The only sensitive parameter is the Oort's constant A, but 
unrealistic values beyond 17 km kpc~^ are required to 
slightly improve the fit. Other scenarii have been studied. 

One possibility was to allow the Belt to cross the 
Galactic plane before reaching its present orientation. The 
effect of an initial rotation has also been considered by 
adding a tangential component to the initial expansion 
velocity. A more realistic braking of the shell has been 
introduced in the late stages of the evolution by grad- 
ually applying a drag force on the shell elements when 
their velocity falls into the 30 to 10 km s~^ interval. In 
parallel, shell fragmentation has been implemented by in- 
creasing the porosity to the interstellar flow. The porosity 
coefficient ap varies linearly from a pure snowplough case 
above 20 km s~^ {ap = 0) to a pure drag without accre- 
tion below 5 km s^^ {ap — 1). The drag force is given 
by dP = —ap Cx Pgas {^' ■ d§) if' where Cx denotes the 
friction coefficient. A value of 0.5 was adopted correspond- 
ing to turbulent friction on a sphere. dS is an elementary 
surface of the shell front, pointing outwards, and if' gives 
the relative velocity of dS with respect to the interstellar 
medium. 

The maximum-likelihood values reached for these dif- 
ferent scenarii are given in Table El Probabilities for ran- 
dom fluctuations from a model to fit the data as well as 
the original expansion model are estimated from the like- 
lihood ratio -2 InA = -2 lii{Cmax model / Cmax exp) be- 
tween the tw o models. This ratio follows a Xi distribution 
l Eadie et al.l 11971 "l. The corresponding probabilities are 
given in Table O It can be seen that a scenario including 
a late shell fragmentation or a Galactic plane crossing can 
hardly reproduce the data. The improvement in the qual- 
ity of the fit of the expansion model over these scenarii is 
very significant. 

Because of the reduced deceleration of the fragmented 
shell, a lower initial velocity is needed to match the low 
cloud velocities and the small size and low inclination re- 
sulting for the evolved Belt in this case explain the much 
poorer fit. 

The crossing of the Galactic plane results in an hour- 
glass like distorsion of the shell because of the interstellar 
density gradient implying different expansion velocities in 
and out of the Galactic disc. A maximum-likelihood solu- 
tion could, however, be found. It requires a larger initial 
inclination to boost the Belt to larger altitudes and make 
use of the increased gravitational pull towards the plane 
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Table 3. Maximum likelihood values reached for different scenarii. P is the probability that random fluctuations from 
a given model fit the data as well as the expansion model. Characteristic Belt parameters are given for a selection of 
scenarii. 





Vr-rif fkm/sl 




P 




Ei flO*'' J) 


/o \ 

^today \ ) 


'-''today 


b+^rlr,,. (DCI 

'-'today Vl-"-"/ 


expansion 





-625.5 




26.4 ±0.4 


1.0 ±0.1 


17.2 ±0.5 


373 ±5 


233 ±5 


rotation 


5 


-626.2 


25.4% 












rotation 


50 


-625.9 


35.9% 












rotation 


-50 


-625.6 


77.7% 


26.4 ±0.5 


1.1 ±0.1 


17.3 ±0.5 


372 ±6 


232 ±6 


rotation 


100 


-626.0 


34.8% 












rotation 


-100 


-625.5 


100% 


26.4 ±0.5 


1.0 ±0.1 


17.1 ±0.5 


371 ±5 


232 ±5 


plane crossing 





-936.3 


< 


51.8 ± 1.0 


0.6 ±0.06 


15.1 ±0.7 


304 ±6 


219 ±6 


fragmentation 





-667.6 


5e-20 


28.9 ±0.6 


0.3 ±0.03 


5.7 ±0.5 


142 ± 5 


72 ±5 



to cross it and reach a 20° tilt on the other side. A lower 
initial kinetic energy is needed because of the extra energy 
provided by the disc torque. The evolved Belt disc ends 
up being quite warped. As expected from Figure the 
resulting fit is much poorer than without plane crossing. 
This is confirmed by a dramatic decrease in the maximum 
likelihood. This scenario requires a substantially longer 
time span of 51.8 ± 1.0 Myr. 

The fit seems rather insensitive to an initial rotation, 
even for large rotation velocities comparable to the ex- 
pansion one. The age and energy estimates, as well as the 
final geometry, are not affected by an initial rotation. This 
set of results show that the dominant factor in the shell 
dynamics and final velocity distribution is the momentum 
accreted from the interstellar gas. 

4.3. Age of the Belt 

The best value found for the age of the Belt, r = 26.4±0.4 
Myr, is comparable to previous estimates b ased on the ga s 
dynamics (31-36 Myr in the 2D model of lOland l|l982l) . 
23 Myr and 1 5.5 M yr in the 3D superbubble model of 
iMoreno et all l)l999|) . without and with interstellar pres- 
sure). The shorter timescales involved in the 3D models 
may be due to the faster expansion in the rarefied medium 
away from the Galactic disc in the early stages. So, the age 
estimate is sensitive to the interstellar stratification. The 
dynamics of the stellar system, mainly its vertical oscilla- 
tions and est imates of the P ort constants, imply an age of 
34 ± 3 Myr ((Comer6nlll999t) . Taken together, these vari- 
ous measurements suggest a dynamical age of the Belt of 
order 30 Myr. 

Yet, stellar ages point to a twice longer timescale for 
the Belt system. Determining the age of the massive stars 
from photometric measuremen ts. iWestin (1985.) p roposed 
an upper limit of 60 Myr, and iTorra et al.l |200o|l showed 
that 60 to 66% of the stars younger than 60 Myr and closer 
than 600 pc belong to the Belt. The X-ray luminosities of 
2q30.5±o.5 g].g^g of the young low- mass stars correspon d 
to ages between 30 and 80 Myr ^Guillout et alJll998l) . 
One difficulty in determining the Belt age from the stars 
is to reliably separate the Belt and the Galactic popula- 
tions. Another comes from deriving the photometric age 
without properly taking into account the stellar rotation. 



iFig^ieras" Blasil l|l998l) indeed found that photometric 
ages of rotating B7-A4 stars are overestimated by 30 to 
50% on average. This important bias would also affect the 
highly rotating massive stars and would help reduce the 
discrepancy between the dynamical age of 30 Myr and 
the stellar age of ~ 60 Myr. On the other hand, an origin 
in a single explosive event or continuous energy injection 
would imply different dynamical ages. As discussed in the 
next section, gradual injection would further reduce the 
dynamical timescale. One way to reduce the dynamical 
and stellar age discrepancy is to allow the Belt to cross 
the Galactic plane before reaching its present orientation. 
The much poorer fit for this scenario, however, does not 
support this possibility. 

4.4. Input energy 

The expansion results from an initial kinetic energy in- 
put of (1.0 ±0.1) lO'*^ J which is comparable, but slightly 
higher than the en ergy of 6.1 1 0^^ J required in the 2D 
expansion model of lOlanol l)l982|) . This is due to the extra 
amount of work needed to expand against the gravita- 
tional pull of the Galactic disc in the early phases. An 
equivalent energy of 6 10*** J is required for the 3D expan- 
sion of a superbubble, in a 3 times lower interstellar den- 
sity, but against ambient pressure (.Moreno et al. 1 9991) . 
The 10 times lower value they found for pressureless ex- 
pansion is due to their choice of a very low density. An 
initial energy of 10^*^ J is typical of the energy deposit 
from multiple supernovae inside a young stellar cluster, 
but it would require a series of explosions within a rather 
short period. The energy budget is also typical of a hy- 
pernova powering a 7-ray burst event. Exploring whether 
the asymetrical explosion can explain the Belt tilt is un- 
der study. The energy budget can also be accounted for 
by the potential energy of high- velocity clouds falling on 
the Galactic disc (jComeron fc Torrall992HComer6n et alJ 
■19941 

Continuous energy injection by multiple supernovae 
and expansion into a highly inhomogeneous medium 
would yield a larger dynamical range in velocity. Injecting 
energy gradually would, however, accelerate the Belt ex- 
pansion. To zcroth order approximation, expansion in a 
uniform medium of density p, with constant power injec- 
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Fig. 7. The Gould Belt evolution as seen at different epochs after the outburst (see Figure [Tj), assuming a crossing of 
the Galactic plane. The plot for 50 Myr is representative of the present-day geometry. 



tion Pinj, can be modelled by -^{npHR'^R) = 2TTHRpint, 
where pint represents an internal source of pressure to 
account for the energy deposit. Its mechanical power is 
2Tr HRpintR = Pinj- Solving for power-law expansion 
(pint oc R~'' and R oc t"), one finds a faster expansion 
R oc f'^/* for constant power injection than from a single 
explosion {R oc t^^^). The faster evolution would further 
reduce the dynamical age with respect to the stellar age, 
and also reduce the Belt eccentricity with respect to the 
location of the clouds and OB associations. The complex 



interplay between the power injection and momentum ac- 
cretion is being investigated. 



4.5. Mass accretion 

The swept- up mass in the evolved shell is found to be 2.4 
lOf_M0. The larger value of 1.2 10*^ M© proposed bv lOland 
ill982^ results from his choice of a 20 % larger density in 
the Galactic plane and the absence of vertical gradients 
and expansion away from the plane. The estimate of 3.3 



PerrotjC.A. and Grenier I.A.: 3D dynamical evolution of the Gould Belt 



11 




150 100 50 -50 -100 -150 

LON 



Fig. 8. Total integrated mass along the Gould Belt rim, 
in solar masses per pc^, in the snowplough case, at five 
different epochs over the past 5 Myr. 

10^ Mq found bv lMoreno et alJ l|l999l) in the 3D supershell 
is close to our result. The 3 times lower density value they 
adopted is compensated by the lack of vertical gradients. 

Figure IHl shows the integrated mass distribution along 
the Belt over the past 5 Myr. The mass accumulated in the 
snowplough case is lowest in the directions 35° < I < 100° 
and —155° < I < —110°. These regions in the first and 
third quadrants are indeed relatively free of major cloud 
complexes (c./. Figure [SJ. Conversely, the swept-up mass 
is highest in the -100° < I < 22° and -165° < Z < 115° 
intervals that encompass most of the main nearby com- 
plexes. The nodal line of maximum accretion is close to 
the Belt nodal line, illustrating the importance of the Belt 
tilt and of the vertical density gradient. In other words, 
the Belt accretes more near the Galactic disc. Figure [S] il- 
lustrates the Belt precession, of the order of a few degrees 
per Myr, that can be attributed to the Galactic differential 
rotation. 

4.6. The Belt wave and the OB associations 

The velocit y field of the nearby y oung stellar groups (see 
Figure 5 in lLindblad et alJ l)l997j) ') suggests a stream-like 
motion spreading out of the Belt region. Because of the 
frictionless motion of the stars in the Galactic potential, 
one expects older stellar groups, born in a faster expand- 
ing Belt shell, to travel farther out than younger groups 
born from a slowed-down Belt. The opposite situation is 
observed in Figure |51 where younger OB associations are 
found at larger radii than older ones (notice in particular 
the LCC/UCL/US and the Oril a/b/c sequences). The 
present-day positions and velocities of stars born 10 Myr 
ago from the expanding shell have been computed from 
their ballistic motion in the Galactic potential. Note that 
velocities were corrected for solar motion (using Uxq — 10 
km s~^, VyQ — 5 km s~^, and Wzq = —7 km s^^) and 
for Galactic rotation (using Oort's constants A = 14.8 km 




-400 -300 -200 -100 100 200 300 

Y(pc) 



Fig. 9. Present-day velocities (projected on the Galactic 
plane) of the nearby OB associations (squares) 
jde Zeeuw et all Il999l) . and of stars that would have 
been born in the expanding Belt shell 10 Myr ago. 
Their position (circles) and velocity today have been 
computed after ballistic motion in the Galactic potential. 
All velocities are corrected for solar motion and Galactic 
rotation. Filled and open symbols mark negative and 
positive vertical Vz velocities, respectively. The presumed 
age for each OB association is given in Myr. The Sun is 
marked by a star. 

s~^ kpc~^ and B = —12.4 km s^^ kpc~^) in order to allow 
a direct comparison with Figure 29 from Ide Zeeuw et 
(1299). It implies that velocities are considered with re- 
spect to the Local Standard of Rest at the position of the 
object considered. One expects the computed present-day 
positions to be systematically too far out because the Belt 
shell model has been fitted against the current position 
of the OB associations, not against their birth positions 
several Myr ago. There is, however, no convincing corre- 
lation between the observed and expected velocity fields, 
particularly in the vertical velocity component due to the 
gravitational pull from the Galactic disc. The current data 
suggests that the Belt wave triggers star formation when 
overtaking a cloud and that the association average veloc- 
ity does not relate to that of the progenitor cloud. 

5. Conclusion 

The Gould Belt expansion has been modelled in 3D and 
has been fitted against the current direction and velocity 
of the H2 and HI clouds that are found in the solar neigh- 
bourhood. Distance measurements have also been used for 
a subset of clouds. Due to the combined effects of the 
Galactic differential rotation, its gravitational pull, and 
the interstellar density gradients, the present Belt section 
has evolved into a slightly warped ellipse, with semi-axes 
of (373 ± 5) pc and (233 ± 5) pc, and an inclination of 
17.2° ± 0.5°. Its centre is currently located (104 ± 4) pc 
away from the Sun, in the direction Icentre — 180.4° ±2.0°. 
The thickness that best fits the data is if = (60 ± 1) 
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pc. While these characteristics nicely compare with previ- 
ous estimates, a different Belt orientation has been found 
that is driven by the presence of new major H2 com- 
plexes and the revised distance information used here from 
Hipparchos measurements. The Belt position and orienta- 
tion arc foimd to coincide with the main OB associations 
and molecular clouds in our vicinity. A larger velocity am- 
plitude would better represent the gas pattern, but it turns 
out to be insensitive to all model ingredients but for the 
Oort constant A. Other scenarii have been tested including 
shell fragmentation in the late stages, an initial rotation, 
or the crossing of the Galactic disc. They did not yield 
an increase in the final dynamical range in velocity. Much 
poorer fits to the data were obtained for the fragmented 
shell and the disc crossing cases. The Belt evolution has 
been compared to the kinematics of the nearby OB as- 
sociations. Its expansion appears to have no influence on 
the observed average stellar motions. On the other hand, 
the swept-up mass distribution along the Belt rim, for 
a total of 2.4 10^ M©, is reasonably consistent with the 
cloud distribution. The younger OB associations are sur- 
prisingly found farther away from the Belt centre than the 
older ones. 

The initial kinetic energy of Ei = (1.0 ± 0.1) lO"*'"^ J 
amounts to that of 10 supernovae, as previously proposed 
by several authors. The dynamical age of (26 ±0.4) Myr is 
equivalent to earlier estimates based on the gas expansion, 
but it is only half of the 60 Myr age derived from photo- 
metric stellar ages. This discrepancy does not result from 
the choice of a single explosion at the origin. Continuous 
energy injection would increase the discrepancy by further 
reducing the dynamical age. The age estimate is not sen- 
sitive to an initial rotation either. Allowing one crossing of 
the Galactic disc before reaching the present inclination 
implies a dynamical age of 52 Myr in better agreement 
with the stellar ages, but the very poor fit docs not sup- 
port this possibility. Important biases in the photometric 
derivation of stellar ages for rapidly rotating stars have 
been reported that could help solve this discrepancy. 

Acknowledgements. The authors would like to thank J. de 
Bruijne and A. Brown for kindly providing data on the OB as- 
sociations. Wo gratefully acknowledge the referee, A. Blaauw, 
for his helpful comments and discussions. We also thank J.-P. 
Chieze for useful discussions on hydrodynamical issues. 

References 

Comeron, F. 1999, A&A, 351, 506 
Comeron, F. & Torra, J. 1992, A&A, 261, 94 
— . 1994, A&A, 281, 35 

Comeron, F., Torra, J., & Gomez, A. E. 1994, A&A, 286, 
789 

Crczc, M., Chcreul, E., Bienayme, O., & Pichon, C. 1998, 

A&A, 329, 920 
Dame, T. M., Ungcrechts, H., Cohen, R. S., et al. 1987, 

ApJ, 322, 706 

de Geus, E. J., Bronfman, L., & Thaddeus, P. 1990, A&A, 
231, 137 



de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J., 

Brown, A. G. A., & Blaauw, A. 1999, AJ, 117, 354 
Dickey, J. M. & Lockman, F. J. 1990, ARA&A, 28, 215 
Eadie, W. T., Drijard, D., James, F. E., Roos, M., & 

Sadoulet, B. 1971, Statistical Methods in Experimental 

Physics (Amsterdam, North-Holland) 
Figueras, F. & Blasi, F. 1998, A&A, 329, 957 
Gehrels, N., Macomb, D. J., Bertsch, D. L., Thompson, 

D. .1., & Hartman, R. C. 2000, Nature, 404, 363 
Gould, B. A. 1874, in Proceedings of the American 

Association for Advanced Science, 115 
Grenier, I. A. 2000, A&A, 364, L93 

Grenier, I. A., Lebrun, F., Arnaud, M., Dame, T. M., & 
Thaddeus, P. 1989, ApJ, 347, 231 

Guillout, P., Sterzik, M. F., Schmitt, J. H. M. M., Motch, 
C, & Neuhauser, R. 1998, A&A, 337, 113 

Hartmann, D. & Burton, W. B. 1997, Atlas of galactic 
neutral hydrogen (Cambridge; New York: Cambridge 
University Press, ISBN 0521471117) 

Herschel, J. F. W. B. 1847, Results of astronomical obser- 
vations made during the years 1834, 5, 6, 7, 8, at the 
Cape of Good Hope (London, Smith, Elder and co.) 

Lindblad, P. O. 1967, Bull. Astron. Inst. Netherlands, 19, 
34+ 

Lindblad, P. O., Palous, J., Loden, K., & Lindegren, L. 
1997, in ESA SP-402: Hipparcos - Venice '97, Vol. 402, 
507-512 

Maddalena, R. J., Moscowitz, J., Thaddeus, P., & Morris, 

M. 1986, ApJ, 303, 375 
Moreno, E., Alfaro, E. J., & Franco, J. . 1999, ApJ, 522, 

276 

Murphy, D. C., Cohen, R., & May, J. 1986, A&A, 167, 

234 

Murphy, D. C. & Myers, P. C. 1985, ApJ, 298, 818 
Ncldcr, J. A. & Mead, R. 1965, Comput. J., 7, 308 
Ojha, D. K., Bienayme, O., Robin, A. C, Creze, M., & 

Mohan, V. 1996, A&A, 311, 456 
Olano, C. A. 1982, A&A, 112, 195 
— . 2001, AJ, 121, 295 

Poppel, W. G. L. & Marronetti, P. 2000, A&A, 358, 299 
Poppel, W. G. L. 1997, Fundamentals of Cosmic Physics, 

18, 1 

Sfeir, D. M., Lallement, R., Crifo, F., & Welsh, B. Y. 1999, 

A&A, 346, 785 
Strong, A. W. 1985, A&A, 150, 273 

Strong, A. W., Riley, P. A., Osborne, J. L., & Murray, 

J. D. 1982, MNRAS, 201, 495 
Taylor, D. K., Dickman, R. L., & Scoville, N. Z. 1987, 

ApJ, 315, 104 

Torra, J., Fernandez, D., & Figueras, F. 2000, A&A, 359, 
82 

Ungerechts, H. & Thaddeus, P. 1987, ApJS, 63, 645 
Westin, T. N. G. 1985, A&AS, 60, 99 
Williams, J. P., de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 
693 



PerrotjC.A. and Grenier I.A.: 3D dynamical evolution of the Gould Belt 



13 



Appendix A: CLUMPFIND results 

For each clump detected with CLUMPFIND, the position 
and velocity of the corresponding centroid are given in 
the following Tables. For known molecular complexes, dis- 
tances that were used in the likelihood function (c./. sec- 
tion are also indicated. 



Table A.l. HI clumps 



1 b vlsr D complex 



(deg) 


(deg) 


(km/s) 


(pc) 




98,5 


8,5 


-1,030529 


? 




106 


7 


-6,18318 


? 




16,5 


13,5 


1,030531 


? 




17,5 


12,5 


1,030531 


? 




102,5 


10 


-7,21371 


? 




-169,5 


-13,5 


4,122121 


340 


Ori 


-159 


-9,5 


7,213711 


500 


Ori^ 


-173 


-14,5 


3,091591 


340 


Ori 


-157,5 


-13,5 


6,183181 


500 


Ori^ 


-140 


-12,5 


4,122121 


? 




-168,5 


-13,5 


4,122121 


340 


Ori 


-161 


-8 


7,213711 


? 




-160,5 


-9 


6,183181 


? 




68 


-8 


8,244241 


? 




-76 


6 


-4,943996 


? 




-77,5 


5 


-3,295996 


? 




-17 


8 


4,120004 


? 




-18,5 


9 


3,296004 


140 


Lup*^ 


-24 


8 


4,944004 


140 


Lup'' 


-1 


8 


4,120004 


125 


Oph'^ 


-21,5 


9 


5,768004 


140 


Lup*" 


-15,5 


8 


4,120004 


? 




-20,5 


8 


2,472004 


? 




-29,5 


8 


4,944004 


? 




-26,5 


9 


4,120004 


? 




-30,5 


8 


4,944004 


? 




11 


9 


5,768004 


? 




-58 


5 


1,648004 


? 




-34,5 


9 


1,648004 


? 




9 


8 


3,296004 


? 




-42 


6 


0,000004 


? 




-93,5 


-9 


3,296004 


380 


Vela 


-97 


-8 


9,888004 


380 


Vela 


-92 


-9 


3,296004 


380 


Vela 


-42,5 


-5 


-4,943996 


7 




-83,5 


-9 


2,472004 


7 




-72 


-6 


-3,295996 


7 




-42,5 


-5 


-2,471996 


7 




-69 


-5 


-2,471996 


7 




-85,5 


-8 


0,824004 


7 




-69 


-6 


-4,119996 


7 




-77 


-10 


0,000004 


7 




-66,5 


-5 


-4,119996 


7 




-61,5 


-5 


-2,471996 


7 




-50 


-8 


-0,823996 


7 




-16,5 


-9 


-2,471996 


7 




-55 


-10 


1,648004 


7 





Maddalena et all \l98dt 


Murohv et al.' 


(■19861) 


de Gcus ct al. 


(,1990.) 
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Table A.2. CO clumps 


1 


b 


^ LSR 


D 


complex 


(deg) 


(dee) 


(km/s) 


(pc) 




11 


-8 


4,1466 


? 




133,5 


9,5 


-4,5462 


350 


Cas^ 


114 


14 


-5,8466 


300 


Cep^' 


-26 


17,5 


4,5566 


? 




-170,5 


-14 


-0,645 


340 


Ori 


101,5 


15 


1,9558 


300 


Cep'' 


-154 


-15 


4,5566 


500 


Ori'' 


108 


16,5 


-4,5462 


300 


Cep^' 


161 


-9 


-4,5462 


320 


Per'^ 


-22 


22 


3,2562 


? 




-19,5 


8,5 


3,2562 


140 


Lup"* 


125 


12 


0,6554 


350 


Gas'' 


5,5 


19,5 


3,2562 


? 




93 


9 


-3,2458 


? 




117 


9 


-3,2458 


? 




133,5 


11 


-3,2458 


350 


Gas'' 


44,5 


8,5 


4,5566 


? 




-155,5 


-12,5 


-0,645 


? 




-143,5 


-16,5 


3,2562 


? 




-154 


-17 


8,4578 


500 


Ori'' 


-150,5 


-20 


7,1574 


500 


Ori'' 


-7,5 


16 


1,9558 


125 


Oph'^ 


-7 


15 


3,2562 


125 


Oph'^ 


-149,5 


-20 


4,5566 


? 




-155,5 


-14,5 


8,4578 


500 


Ori'' 


-21,5 


15,5 


3,2562 


? 




158 


-21 


5,857 


320 


Per'^ 


-148 


-19,5 


1,9558 


? 




-18 


9 


3,2562 


140 


Lup"* 


-5 


20 


1,9558 


125 


Oph'= 


-4 


16 


0,6554 


125 


Oph" 


-22,5 


16 


4,5566 


? 




-5,5 


14 


3,2562 


125 


Oph'^ 


-168,5 


-11,5 


8,4578 


? 




-21,5 


14,5 


3,2562 


? 




159 


-20,5 


5,857 


320 


Per'' 


-4 5 


18,5 


-0,645 


? 




160 


-18,5 


7,1574 


? 




-5,5 


15,5 


0,6554 


125 


Oph'= 


-3 


18,5 


0,6554 


125 


Oph'= 


-22 


17 


4,5566 


? 




-165,5 


-17 


-3,2458 


340 


Ori 


-1,5 


21 


-0,645 


? 




0,5 


9,5 


3,2562 


125 


Oph'= 


112,5 


16 


-5,8466 


300 


Gep'' 


-1,5 


18 


0,6554 


125 


Oph'= 


1 


8,5 


1,9558 


? 




165 


-9,5 


-1,9454 


140 


Tau^ 


103,5 


13,5 


0,6554 


300 


Gep" 


-154,5 


-8,5 


8,4578 


500 


Ori'' 


-143 


-14,5 


8,4578 


? 




20,5 


9,5 


5,857 


250 


Aqu 



a 


Grenier et alJ 


(119881) 


b 


Maddalena et al." (l98&) 


c 


Uneerechts & ThaddeuJ 


d 


Murohv et al. 




e 


de Geus et al. 


I1990I! 


f 


Murohv & Mverd (Il985l) 



Table A. 3. CO clumps (continued) 



1 


b 


^LSR 


D 


complex 


(dee) 


(deg) 


(km/s) 


(pc) 




132,5 


7 


-4,5462 


350 


Gas'' 


-158,5 


-11,5 


7,1574 


500 


Ori'' 


158,5 


-9 


-5,8466 


320 


Per'' 


133,5 


9,5 


-4,5462 


350 


Gas 


114 


14 


-5,8466 


300 


Gep'' 


-26 


17,5 


4,5566 


? 




-170,5 


-14 


-0,645 


340 


Ori 


101,5 


15 


1,9558 


300 


Gep'' 


-154 


-15 


4,5566 


500 


Ori" 


108 


16,5 


-4,5462 


300 


Gep'' 


161 


-9 


-4,5462 


320 


Per'' 


-22 


22 


3,2562 


? 




-19,5 


8,5 


3,2562 


140 


Lup" 


125 


12 


0,6554 


350 


Gas 


5,5 


19,5 


3,2562 


? 




93 


9 


-3,2458 


? 




117 


9 


-3,2458 


? 




133,5 


11 


-3,2458 


350 


Gas 


44,5 


8,5 


4,5566 


? 




-155,5 


-12,5 


-0,645 


? 




-143,5 


-16,5 


3,2562 


? 





